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ON THE BACKWARD STABILITY OF THE 
SCHWARZSCHILD BLACK HOLE SINGULARITY 


GRIGORIOS FOURNODAVLOS 


Abstract. We study the backwards-in-time stability of the Schwarzschild sin¬ 
gularity from a dynamical PDE point of view. More precisely, considering a 
spacelike hypersurface Sq in the interior of the black hole region, tangent to 
the singular hypersurface {r = 0} at a single sphere, we study the problem 
of perturbing the Schwarzschild data on Sq and solving the Einstein vac¬ 
uum equations backwards in time. We obtain a local well-posedness result for 
small perturbations lying in certain weighted Sobolev spaces. No symmetry 
assumptions are imposed. The perturbed spacetimes all have a singularity at 
a “collapsed” sphere on Sq, where the leading asymptotics of the curvature 
and the metric match those of their Schwarzschild counterparts to a suitably 
high order. As in the Schwarzschild backward evolution, the pinched initial 
hypersurface Eq ‘opens up’ instantly, becoming a smooth spacelike (cylindri¬ 
cal) hypersurface. This result thus yields classes of examples of non-symmetric 
vacuum spacetimes, evolving forward-in-time from smooth initial data, which 
form a Schwarzschild type singularity at a collapsed sphere. We rely on a 
precise asymptotic analysis of the Schwarzschild geometry near the singularity 
which turns out to be at the threshold that our energy methods can handle. 
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1. Introduction 

It is well-known (cf. Birkhoff’s theorem [10) 1 that the only spherically symmetric 
solution to the Einstein vacuum equations (EVE) 

(1.1) R.iCa&(ff) = 0, 

is the celebrated Schwarzschild spacetime. It was in fact the first non-trivial so¬ 
lution to the Einstein field equations to be discovered [10]. In Kruskal (null) u, v 
coordinates the maximally extended metric reads 

(1.2) = —U^dudv + r'^{d9^ + sin^ Odcj)^), 

where 12^ = e~™, M > 0, and the radius function r is given implicitly by 

r T 

(1.3) uv = il - )e™. 

Here the underlying manifold is endowed with the differential structure of 



Figure 1. The Kruskal plane. 
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U X S^, where U is the open subset {uv < 1} in the uv plane; see Figure The 
spacetime has an essential curvature singularity at r = 0, (the future component of) 
which is contained in the interior of the black hole region, the quadrant u > 0, ?; > 0. 
In fact, a short computation shows that the Gauss curvature of the rtw-plane equals 


(1.4) 


■^K = 


2M 

j-3 


and hence the manifold is inextendible past r = 0. An interesting feature of 
this singularity is its spacelike character, that is, it can be viewed as a spacelike 
hypersurface. 


Yet another interesting feature of the Schwarzschild singularity is its unstable nature 
from the evolutionary dynamical point of view. To illustrate this consider a global 
spacelike Cauchy hypersurface in Schwarzschild (Figure]^. An initial data set 
for the EVE consists of a Riemannian metric g on S and a symmetric two tensor 
K verifying the constraint equations 


Kij — WitlgK = 0 
R-\K\^ + {tTgK)^ =0 ’ 


where V,R are the covariant derivative and scalar curvature intrinsic to g. 



Figure 2. 


The instability of the Schwarzschild singularity (w.r.t. the forward Cauchy prob¬ 
lem) can already be seen by examining the maximal developments of initial data 
sets on S arising from the celebrated Kerr [5] (explicit) 2-parameter JC{a,M) fam¬ 
ily of solutions - of which Schwarzschild is a subfamily (a = 0). For a yf 0 the 
singularity completely disappears and the corresponding (maximal) developments 
extend smoothly up to (and including) the Cauchy horizons (which lie far beyond 
the hypersurface uv = 1). Moreover, taking |a| <C 1, the ‘difference’ of the cor¬ 
responding initial data sets from the Schwarzschild one (with the same M > 0), 
measured in standard Sobolev normsQcan be made arbitrarily small. 


^The difference can be defined, for example, component wise for the two pairs of 2-tensors 
with respect to a common coordinate system and measured in Sobolev spaces used in the 

literature [4]. 
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In fact, the Schwarzschild singularity is conjecturally unstable under generic pertur¬ 
bations on S. According to a scenario proposed by Belinskii-Khalatnikov-Lifshitz 
[2] originally formulated for cosmological singularities, in general, one should ex¬ 
pect solutions to exhibit oscillatory behaviour towards the singularity^ Although 
this work was not supported by much rigorous evidence, it has received a lot of 
attention over the years, see [25] and the references therein (and [T^ for related 
numerics). On the other hand, there is a growing expectation that, at least in a 
neighbourhood of subextremal Kerr, the dominant scenario inside the black hole 
is the formation of Cauchy horizons and (weak) null singularities. This has been 
supported by rigorous studies on spherically symmetric charged matter models, see 
works by Poisson and Israel |2l], Ori [23| and recently by Dafermos |9|. 

However, it is not inferred from the existing literature whether the non-oscillatory 
type of singularity observed in Scwarzschild is an isolated phenomenon for the EVE 
in some neighbourhood of the Schwarzschild initial data on E or part of a larger 
family. A priori it is not clear what to expect, since one might argue that such a 
special singularity is a mathematical artefact due to spherical symmetry. Therefore, 
we pose the following question: 

Is there a class of non-spherically symmetric Einstein vacuum spacetimes which 
develop a first singularity of Schwarzschild type ? 


The goal of the present paper is to answer the preceding question in the affirma¬ 
tive. A Schwarzschild type singularity here has the meaning of a first singularity 
in the vacuum development which has the same geometric blow up profile with 
Schwarzschild and which can be seen by a foliation of uniformly spacelike hyper¬ 
surfaces; hence, not contained in a Cauchy horizon. We confine the question to the 
formation of one singular sphere in the vacuum development in the same manner 
as in Schwarzschild, where each point on the sphere can be understood as a distinct 
ideal singular point of the spacetime in the language of TIPs [12] . Ideally one would 
like to study the forward problem and identify initial data for the EVE on E (Figure 
that lead to such singularities. Although this is a very interesting problem, we 
find it far beyond reach at the moment. Instead we study the existence problem 
backwards-in-time. 


More precisely, we adopt the following plan: Let Eg be a spacelike hypersurface in 


at a single sphere of the singular hyperbola r = 0 inside 
We assume, without loss of generality!^ that the tangent 


Schwarzschild, tangem 
the black hole; Figure [3 
sphere is (u = 1, u = 1) in Kruskal coordinates (1.21. Consider now initial data sets 
{g, K) on Eg for the EVE (1.11 , which have the same singular behaviour to leading 
order at (u = l,u = 1) with the induced Schwarzschild initial data set {^g,^K) 
on Eg and solve the EVE backwards, as depicted in the 2-dim Figure [^ without 
symmetry assumptions. 


^Homogeneous spacetimes with singularities that exhibit oscillatory behaviour have been rig¬ 
orously studied by Ringstrdm m for the Euler-Einstein system with Bianchi IX symmetry. 

^The tangency here should be understood with respect to the differ entia l structure of the 
Kruskal maximal extension induced by the standard u, v, 6, 0 coordinates ( |l.2[ |. 

^Recall that the vector field tangent to the r = const, hypersurfaces (Figure Q is Killing and 
we may hence utilize it to shift Eq and (u = l,u = 1) to whichever point on {uv = 1} we wish; 
Figure (j^. 
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Realizing the above plan we thus prove the existence of a class of non-spherically 
symmetric vacuum spacetimes for which (1) the leading asymptotics of the blow up 
of curvature and in general of all the geometric quantities (metric, second funda¬ 
mental form etc.) coincide with their Schwarzschild counterparts, as one approaches 
the singularity, and (2) the singularity is realized as the limit of uniformly space¬ 
like hypersurfaces, which in the forward direction “pinch off’ in finite time at one 
sphere. Conversely, we visualise the backward evolution of {Tiq, g,K) in the follow¬ 
ing manner: At ‘time’ t = 0 the initial slice Eq is a two ended spacelike (3-dim) 
hypersurface with a sphere singularity at (it = 1, u = 1). Once Eg evolves through 
(1.1), it becomes instantaneously a smooth spacelike hypersurface ET-,r > 0 and 
the singular pinch opens up; Figure]^ 



Figure 4. 


The main difficulty to overcome in the backward local existence problem is the 
essential singularity on Eg, which of course renders it beyond the scope of the 
classical local existence theorem for the Einstein equations [3], even its latest state of 
the art improvement by Klainerman-Rodnianski-Szeftel m, which requires at the 
very least the curvature of the initial hypersurface to be in . For the Schwarzschild 
initial data set {^g, ^K) on Eg, and hence for perturbed initial data sets (g, K) with 
the same leading order geometry at (it = 1, ?; = 1), it is not hard to check (fj^ that 
the initial curvature is at the singular level 



( 1 . 6 ) 


R^LP(Eo), ViF^LP(Eg) 
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To our knowledge, the only general local existence results, without symmetry as¬ 
sumptions, for the EVE (1.1) with singular initial curvature not in have been 


achieved only fairly recently by Luk-Rodnianski HZIIIH] and Luk m for the char¬ 
acteristic initial value problem. However, the singular level of the Schwarzschild 
singularity is much stronger than the ones they consider; delta curvature singu¬ 
larities [mils] and weak null singularities [16) . Also, since we are interested in 
the spacelike character of the singularity of the perturbed spacetimes we wish to 
construct, the non-characteristic initial value problem seems to be a more natural 
setting to consider. 

We proceed now to formulate a first version of our main results; for more precise 


statements, in terms of weighted Sobolev spaces, see Theorems 4.7 4.9 6.3 


Theorem 1.1. There exists a > 0 sujficiently large, such that for every triplet 
{^o, 9 ,K) verifying: 


’AT -|- r“u, where 0,u are 2-tensors on Eq bounded in 


(i) the constraints (1.5), 

\ii) g = ^g r°‘0, K = 
respectively, 

{Hi) ||5- ^5lU“(So) < 1, 

there exists a H'^ local solution g to the Einstein vacuum equations (1.1) with ini¬ 


tial data {g,K), unique up to isometry, in the backward region to Sq, foliated by 
{ET-}T-g[o,T] (Figure^; the time of existence T > 0 depends continuously on the 
norms of 0,u and the exponent a > 0. 


The fact that non-trivial initial data sets in compliance with Theorem o exist is 
not at all obvious nor standard. It has to be shown essentially that for any large 
parameter a > 0, the set of non-spherically symmetric solutions to the constraint 
equations (1.5), having the asymptotics {ii), is non-empty. This is achieved in the 
following theorem. 


Theorem 1.2. Let a > 0 be sufficiently large, consistent with Theorem \l . 1\ There 
exists a small ball B C H'^{d'Eo) x H^{dT,Q), such that for every pair {0,u) of 
2-tensors defined on the boundary of Eg, which lie inside that ball, there exists 
at least one extension of {O, u) as 2-tensors in Eg satisfying the assumptions of 
Theorem \l.l\ In particular, the pair (g,K) defined via {ii) in Theorem \l.l\ satisfies 
the constraints ljl.5^. 


Finally, let us emphasize the fact that the above spacetimes are very special in that 
they agree with Schwarzschild at the singularity to a high (but finite) order - this 
is captured by the large exponent a > 0 in Theorem o- and therefore are non¬ 
generic. The need to choose a large may be seen however natural to some extend 
in view of the instability of the Schwarzschild singularity; from the point of view of 
the forwards-in-time problem. Indeed, the stable perturbations of the Schwarzschild 
singularity must form a strict subclass of all perturbed vacuum developments. 


1.1. Method of proof and outline. The largest part of the paper is concerned 
with the evolutionary part of the problem, i.e., proving Theorem |l.l[ Due to the sin¬ 
gular nature of backward existence problem described above, Figure]^ the choice of 
framework must be carefully considered. The standard wave coordinates approach 
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[3] does not seem to be feasible in our situation; one expects that coordinates would 
be highly degenerate at the singularity. Also, the widely used CMC gauge condi¬ 
tion is not applicable, since the mean curvature of the initial hypersurface Eg blows 
up (^. Instead, we find it more suitable to use orthonormal frames and rewrite 
the EVE one order higher as a quasilinear Yang-Mills hyperbolic system of equa¬ 
tions [211 [n], under a Lorentz gauge conditionlj for the corresponding connection 
1-forms. We recall briefly this framework in © 


However, even after expressing the EVE in the above framework, the singular level 
of initial configurations do not permit a direct energy estimate approach. In ad¬ 
dition to (1.6), one can see (©) that neither is the second fundamental form in 


(1.7) 




Note that the latter is at the level of one derivative in the metric. Hence, near 
the singularity the perturbed spacetimes we wish to construct do not even make 


sense as weak solutions of the EVE (1.1). Therefore, it is crucial that we use the 


background Schwarzschild spacetime to recast the evolution equations in a new 
form having more regular initial data. We do this in by considering a new 
system of equations for the ‘difference’ between the putative perturbed spacetime 
and Schwarzschild. The resulting equations have now regular initial data and they 
are eligible for an energy method, but there is a price to pay. The coefficients of 
the new system will depend on the Schwarzschild geometry and will necessarily 
be highly singular at r = 0. We compute in ^ the precise blow up orders of 
the Schwarzschild connection coefficients, curvature etc. Nevertheless, the issue of 
evolving singular initial data has become the more tractable problem of finding 
appropriate weighted solution spaces for the final singular equations. 


In §4.2| we introduce the weighted Sobolev spaces which yield the desired flexi¬ 
bility in proving energy estimates. The right weights are given naturally by the 
singularities in the coefficients of the resulting equations, namely, powers of the 
Schwarzschild radius function r with a certain analogy corresponding to the order 
of each term. After stating the general local existence theorems in §4.3| and a more 
precise version of Theorem o we proceed to its proof via a contraction mapping 
argument which occupies Section Therein we derive the main weighted energy 
estimates by exploiting the asymptotic analysis at r = 0 of the Schwarzschild com¬ 
ponents (©). It is necessary in our result that the power of r, a > 0, in the weighted 
norms is sufficiently large; cf. assumption {ii) in Theorem o In the estimating 
process certain critical terms are inevitably generated, because of the singularities 
in the coefficients of the system we are working with; these terms are critical in 
that they appear with larger weights than the ones in the energy we are trying to 
control and thus prevent the estimates from closing. The exponent a > 0 is then 
picked sufficiently large such that these critical terms have an overall favourable 
sign; this allows us to drop the critical terms and close the estimates. 


The largeness of a forces the perturbed spacetime to agree asymptotically with 
Schwarzschild to a high order at the singularity. Although the latter may seem 
restrictive, it is quite surprising to us that there even exists a suitable choice of 


^The analogue of a wave gauge for orthonormal frames. 
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a which makes the argument work in the first place. A closer inspection of our 
method reveals that it is very sensitive with respect to certain asymptotics of the 
coefficients in the equations that happen to be just borderline to allow an energy- 
based argument to close. The most important of these are the blow up order of 


the sectional curvature (1.4) and the rate of growth of the Schwarzschild radius 
function r backwards in time. The latter corresponds to the ‘opening up’ rate of 
the neck pinch of the singular initial hypersurface Sqj Figure]^ In this sense the 
Schwarzschild singularity is exactly at the threshold that our energy-based method 
can tolerate. 


In the last section, ^ we study the constraint equations ( |1.5[ ) in a perturbative 
manner about the Schwarzschild singular initial data set {^g, ^K) on Eg. First we 
prescribe boundary datcH at cIEq close to Schwarzschild (in usual Sobolev norms, 


see Theorem 1.2) and then look for solutions (g, AT) to constraints inside the class 
of singular initial data sets that are compatible with our assumptions on the evo¬ 
lutionary problem; cf. Theorem |1.1| It is very likely that the generic behaviour 
at (u = 1,?; = 1) of solutions to ( |1.5[ ) with this perturbed data at the boundary 
of Eg will not agree with Schwarzschild to a high enough order for Theorem o 
to be applicable. However, we prove the existence of solutions for which the as¬ 
ymptotic behaviour required by Theorem |l.l| indeed holds. We obtain Theorem |1.2| 
as an application of the implicit function theorem. The proof is based on proper¬ 
ties of the constraint map [20], defined in the weighted Sobolev spaces introduced 
above. In order to prove the surjectivity of the linearized constraint map around 
Schwarzschild, we use an indirect argument from |8] adapted to our context. One 
of the main ingredients in the proof involves the Fredholm property of a linearized 
operator, see Proposition |6.6[ which we prove in §6.2| by deriving certain weighted 
elliptic estimates. 


1.2. Final Comments; Possible applications. To our surprise the present evo¬ 
lution bears some resemblance at an analytical level with a prior work on the stabil¬ 
ity of singular Ricci solitons [Tj. Although of different nature, hyperbolic/parabolic 
(respectively), they share a couple of key features such as the opening up rate of 
the singularity and the “borderline” singularities in the coefficients involved. 

The understanding of the question of stability of singularities in Einstein’s equa¬ 
tions and the behaviour of solutions near them is of great significance in the field. 
However, in general very little is known. In terms of rigorous results, substantial 
progress has been made under symmetry assumptions in the presence of matter |6] 
mmm- Moreover, certain matter models enjoy the presence of a monotonic 
quantity, which has been employed to study the stability of singularity formation 
in the general non-symmetric regime, cf. recent work of Rodnianski-Speck m on 
the FLRW big bang singularity. This is in contrast with the vacuum case and the 
unstable nature of the Schwarzschild singularity. We emphasize again the fact that 
the method developed herein does not impose any symmetry assumptions nor it 
relies on any monotonicity. It should be noted as well that it does not depend 
on whether the particular singularity type is generic or not. On the contrary, we 
hope that it can serve as a robust new method which can be generally employed 


®The boundary of Eg consists of two components diffeomorphic to S^. 
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to produce classes of examples of non-symmetric singular solutions to the Einstein 
field equations, which until now are only known to exist under special symmetry 
assumptions and for which the general stability question may be out of reach. 

After our treatment of singular initial data containing a single sphere of {uv = 1}, a 
reasonable next step would be to study whether the construction of non-spherically 
symmetric vacuum spacetimes containing an arc of the singular hyperbola (Figure 
is possible or even the whole singularity r = 0. Certainly this is a more restrictive 
question and at first glance not so obvious how to formulate it as a backward initial 
value problem problem for the EVE. However, we hope that the method developed 
herein could help approach this direction. 

Lastly, one could try to perform a global instead of a local construction by con¬ 
sidering a Cauchy hypersurface Eg extending to spacelike infinity. We expect this 
follows readily from the work here, but we do not pursue it further. Perhaps a 
gluing construction could also be achieved. 


1.3. Acknowledgements. I would like to thank my advisor Spyros Alexakis for 
his valuable suggestions and interest in this work. I am grateful to Mihalis Dafermos 
for helpful comments on a preliminary version of the present paper. Thanks also go 
to Yannis Angelopoulos, Jesse Gell-Redman, Niky Kamran, Volker Schlue, Arick 
Shao for useful discussions. The author is supported by an Onassis Foundation 
Fellowship. 


2. The Einstein equations as a quasilinear Yang-Mills system 


The Einstein vacuum equations (1.11, by virtue of the second Bianchi identity, imply 
the vanishing of the divergence of the Riemann curvature tensor. Decomposing 
the latter with respect to an orthonormal frame, which satisfies a suitable gauge 
condition, it results to a quasilinear second order hyperbolic system of equations 
for the connection 1 -forms corresponding to that frame, which bears resemblance 
to the semilinear Yang-Mills [5T]. Recently this formulation of the EVE played 
a key role in the resolution of the bounded curvature conjecture m- In this 
section we express the EVE o in the above setting, which we are going to use to 
directly solve the Cauchy problem. This necessitates some technical details which 
are carried out in Appendix]^ Also, to avoid additional computations we write all 
equations directly in scalar non-tensorial formj^ 


2.1. Cartan formalism. Let be a Lorenzian manifold and let {eo,ei, 

62 ; 63 } be an orthonormal frame; gab ■= 'niab = dia.g(~li Ij Ij !)■ Assume also that 
has the differential structure of S x [0,T], where each leaf E x {r} =: Et- 
is a 3-dim spacelike hypersurface. We denote the connection 1-forms associated to 
the preceding frame by 

(2.1) {Ax)ij '■= g{'^xeijCj) = -{Ax)ji, 


^It will be clear though which are the covariant expressions; see also m- 




10 


GRIGORIOS FOURNODAVLOS 


where V is the g-compatible connection of . Recall the definition of the 

Riemann curvature tensor 


(2.2) R^uij ■— ^ ^j') ■ 

By the former definition of connection 1-forms, using rriab to raise and lower indices, 
we write 

^Ca^b — {-^a'}b b-k- 


Hence, we have 

^ = Ve,,(Ve„ei) - {Ak)i''Be 

= B/^(A^)i Cfc -l- (A^)i (Afj,)k Bd — (Afj,)i, (Ak)i‘^Bc 

Therefore, we get the following expression for the components of the Riemann 
curvature 

(2.3) Rfiuij = Bfj_{A,j)ij — B^{Afj_)ij + {Aij)i {Afj^)kj — {Afj_)i {A^)kj 

~ iAfi)i, {Ak)ij + {A,y)fj^ {Ak)ij 

or setting 

(2.4) {[Af^jA^ij = {A^)i {Ai,)kj — {A,y)i {Af^)kj 
we rewrite 

(2.5) •— Rfi^ij — Bpii^Aj^^ij Bjji^A^^j ([H^,Ht/])2j (H^^),yj (^A^'jij^ 

where by standard convention 

{A[ij)v\ {Ak)ij ■= iA^)i, {Ak)ij — (^!/)/j {Ak)ij. 

In the same manner we compute the covariant derivative of the Riemann tensor: 

aRfivij = Ba{F^y)ij — {Aa)J^ {Fkv)ij — {A„) {F^k)ij 
- {A„),^{F^,)kj - (A,)/(F^,),fe 

( 2 . 6 ) = BcriF^„)ij - (A^)'" ^^{F^]k)ij - {[Act, F^^jij 


Recall the transformation law of the above quantities under change of frames: Let 
{cjIq be an orthonormal frame on such that 

(2.7) e„ = O^Bk 

and let {Ax)ij '■= giyxBi,Bj) be the corresponding connection 1-forms. Then 

(2.8) {Ax),, = OlO^^{Ax)bc + A(0(')0=mb,. 

In addition, from ( |2.7| ) we have 

VxBa=X{Ol)Bk + 0’yxBk 
{Ax)a Bd = X{ 0 '^) 6 k + 0 '^{Ax)k‘^Bd 

or 


( 2 . 9 ) 


X{0[) = iAx)j0^,-0’^{Ax)k'. 
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2.2. V X Ric = 0. Now we proceed by assuming that the curl of the Ricci tensor 
of the metric g vanishes: 

( 2 . 10 ) = 0 , 

where Rab ■= R^iab^■ A direct implication of the (contracted) second Bianchi 
identity is that the divergence of the Riemann curvature tensor satisfies 


( 2 . 11 ) = = 

Thus, it follows from ( |2.6[ ) that 

(2.12) - ([A^F^.])„- = 0 


or by 

(2.13) □ {A^)ij - e^e„{A^,)ij - e^"{[A^,A^])ij - {Ak)ij) 

= (A'^)'=[^(F.]fc)^, + ([A^F^.])y, 


where □ := — e§ + ef + + e§ is the non-covariant box with respect to the frame 

Ci- Since 


^ - (A[^),y] 


(2.13) takes the equivalent form 

^{A^)ij — ei,e^{A^)ij = ek{Afj)ij -I-e^([j4^, ^1^])^- -I-e^((A[^),^] {Ak)ij) 

(2.14) + 

ly, i,j = 0,1, 2, 3. We remark that (|2.14 1 is an equation of scalar functions. 


2.3. Choice of gauge. Note that the preceding equation is not of hyperbolic type. 
We convert (2.14) into a quasilinear hyperbolic system of equations by imposing a 
Lorentz gauge condition on the orthonormal frame {eilpl^ 

(2.15) A = {divAx)ij V^{Af^)ij — {Ax^i^f.^)ij = e^{A^)ij — {A^)^ (Ak)ij, 

where by A^ we denote some quadratic expression in the connection coefficients 
{A^)ij varying in ij. This a freedom one has in choosing the frame e^; see Lemma 
A.l Under (|2.15[), the equation (2.14) becomes the quasilinear second order 


n](A,),, = (Al'^),] ek{A^),j + e>^{[A^,A,]),, + e>^{{A[^),]\Ak),,) 

(2.16) -I- (A^)^j^(Fj,]fc)y -I- ([A'^, Ff^„])ij + e,y(A^) -I- e^({A^)^'‘{Ak)i 


wave type gauge essentially for . The Coulomb gauge is another alternative which is used 
in m- We do not employ it here. 
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2.4. The reduced equations; Initial data for EVE. Following ( 2.11)-(2.16) 

we actually see that the equation 

(2.17) - VjR^i+e^(divAx - 

=: H^ij = (LHS of pl6| ) - (RHS of ( |2l6| ), 

holds true for every Lorentzian metric g and orthonormal frame without 

any additional assumptions or gauge condition. We call Hi,ij = 0, i.e., the system 


(2.16), the reduced equations. We note that even after the gauge fixing, the reduced 


equations are not equivalent to the EVE (1.1), but only imply the vanishing of the 


curl of the Ricci tensor (2.10). However, one may suitably prescribe initial data for 


(2.16) such that they lead to solutions of the EVE and which are consistent with 


the Lorentz gauge condition (2.15). 


Now we address the initial value problem for the reduced equations H^ij = 0 


aiming to the EVE. To solve the equation (2.16) one needs an equation relating 
the evolution of the orthormal frame {ci}^ to that of the connection 1-forms. Let 
do, di,d 2 , do be a reference fram^in Ex [0, T] (do transversal direction). We express 
Ci in terms of da'- 


(2.18) 


e, = Otda 


By virtue of the diffeomorphism invariance of the EVE, we may assume that the 
timelike unit vector of the orthonormal frame {e^lo of the spacetime we solve for 
is Co = do. Doing so we deduce 

do(0“) = Ceo{da{ei)) = Cdoida)ei + da{[do,ei]) = 0\Cd„{da)db + da([eo,ei]), 

where C denotes the Lie derivative and da is the 1-form dual to da- Setting 
[do, db] =: Fjp^jdc we rewrite 

(2.19) do(0“) = -O.'Ffo,] + {A^o).]'Ol- 

Now we proceed to formulate the necessary and sufficient conditions on the initial 


data set of the reduced equations (2.16), coupled to (2.19), such that the corre¬ 


sponding solution yields a solution to the EVE. The following proposition is proved 

inM 


Proposition 2.1. Let {A,A)ij,01 be a solution of {2-16), {2-19), arising from initial 
configurations subject to 

(2.20) {A„),j{T = 0) = -{A,,)ji{T = 0) do{A„)ij{T = 0) = -do(^j/)ji(r = 0) 

Oo“(r = 0) = /o“ 

and 


( 2 . 21 ) 


{divAx)ij — A'^ = 0 




Ricab{g) = 0 ^ e'^(H,),^-e,(H'^),^-([H'^,H,]),^-(Hl'^),] {Ak)^^ = 0 

on Eq. Then the latter solution corresponds to an Einstein vacuum spacetime 
{A4^'^^,g) and furthermore the frame {ei\{^ (2-18) is g-orthonormal, eg = do, and 
satisfies the Lorentz gauge condition {2-15). 


®Not orthonormal or coordinates, simply a basis frame. 
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Remark 2.2. Note that the second part of (2.21) includes the constraints (1.5); 
Rot, = Roo~ = 0, ^ = 1, 2, 3, on Eg. Once we have chosen the orthonormal frame 
initially and the initial data components {Ao)ij(T = 0), which correspond to the 
do derivative of {etlg, then the rest of the initial data set of (2.16) is fixed by the 
condition ( 2.21|), i.e., the Lorentz gauge and the EVE on the initial hypersurface 
Eq; see Lemma [A. 1[ Remark | A.2 1 


3. The Schwarzschild components 

We fix an explicit Schwarzschild orthonormal reference frame and compute the cor¬ 
responding connection coefficients, which we then use to find the leading asymp¬ 
totics of the second fundamental form and curvature of the initial singular hyper¬ 
surface Eq in Schwarzschild. Knowing the precise leading blow up behaviour of 
these quantities is crucial for the study of local well-posedness in the next section. 
Eor distinction, we denote Schwarzschild components with an upper left script ^. 

Let us consider a specific foliation of spacelike hypersurfaces E,-, r € [0,r], for 
the backward problem in a neighbourhood of {u = l,v = 1); Figure]^ For conve- 
nienc let 

(3.1) T.r ■■ -]^{u + v) + l = T (u,?;) e (1 - e, 1-I-e)^, r e [0,T]. 



In temporal and spatial coordinates r, x 
(3.2) dr-.= -du-dy, 


3 ''=■ 3 — 3 

'-'.T ■ '-'U 


X= -{u-v), 


the metric (1.2) takes the form 

(3.3) ^g = + il^dx^ + r^{dd^ -I- sin^ 0d(j)^), 


= 


32M3 


is easy to see that the following leading asymptotics we derive are independent of the 
particular choice of foliation. 
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By (1.3),(3.2) r is related to t,x via 


(3.4) 




2M' 

from which one can derive the following formulas: 
(3.5) 




a Si" 

= 4m"= 


Remark 3.1. The above first two identities yield the leading asymptotics: 


(3.6) 




as T,x 


0 . 


Directly from the form of the induced metric on S,-, 

(3.7) ^g = n^dx^ + r'^{d0^ + sin^ 9dcj)^), 
we compute the corresponding induced volume form 

(3.8) dfj,Sg = sin 9dxdOdcj) = + 0(r^)] sm9dxd9d4> 

and its rate of change along dr using (|3.5|: 

(3.9) 


12M^ 1 

drdgsg = [ (1 - r) + 0(-)]d/iSg. 

Normalizing, we define the Schwarzschild orthonormal frame 


(3.10) 


5n = ^ 


1 d 


n dr 


^ 


1 d 


Q dx 


d2 = 


Id^ 

r M 


ds = 


r sin 9 d(f) 

and the relative connection coefficients ^{Ag)ij = ^g{^\7g^di,dj) associated to it. 
A tedious computatioip] shows that the non-zero components read 




(3,11) 


’(^2)02 — ‘^(^3)03 — 


n 1 — T 
AM r 

il X 


’(^2)12 = *(^3)13 = —77- 


AM r 


’(^3)23 = 


cot 9 


Recall the (spacetime) divergence formula of the connection 1-forms X —>■ ^{Ax)ij 
(3.12) ^(divAx)., := d>^^{A,),g - ^Ax.^o^h = - ^(A'^)/^(A,)„ 


^One may calculate the connection coefficients using the Koszul formula 




^g{[dti,di],dj) - ^gi[d^,dj],dg) + ^g{[dj,dg],di) 
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Utilizing (3.51 and (3.11), we check that the first order term in the RHS of (3.12) 
vanishes 


(3.13) 
leaving 

(3.14) 


a'^^(A^),, = 0, 


’{divAx)ij — '®(^3)23’^(^2)ij- 


Remark 3.2. Thus, the orthonormal frame (3.10) satisfies a Lorentz gauge type 
condition (2.15). 

Remark 3.3. Summarizing the above identities and formulas we obtain the fol¬ 
lowing leading asymptotics at r = 0 : 

1 1 „ „ 1 


do 




-r^d-r 


di 


4V2M 


(3.15) 






3-r^dx 

2 

C 

..(Ai+l) ^ ’ 


where C depends on M > 0 and k. Notice that the latter asymptotics are sharp for 
k = 0 and when d^^^ = dg'^\ In fact, the components of the second fundamental 
form of the slices ^Ka = ^{Ai)oi, i = 1,2,3, are exactly at this level. In more 
geometric terms we have (up to constants) 


(3.16) 


T 2 


Itrs/itrl- 3 


’R| ~ 4. 


Thus, employing ( |3.8| ),( |3.6[ ) for r = 0, we see that both the scalar curvature and 
the second fundamental of the initial singular hypersurface Eq are far from being 
square integrable 


’it'pd/rs 


-xX^dx= 


1 


-dx = -boo 


(3.17) 


'o a ;2 
1 


’Rfd/isg ~ / -xx^dx= / — da; =-boo 


/o X 2 


The same holds for the mean curvature of Eg. In fact, a similar calculation shows 
tvs/K ^ LP,p> f. 


4. The local-in time well-posedness 


4.1. Perturbed spacetime; A transformed system. Let {g,K) be a perturba¬ 
tion of the Schwarzschild initial data set {^g,^K) on Eg, verifying the constraints 
(|l.5|), and let \ei}^ be an orthonormal frame of (Eg, 5 ). We fix a reference frame 
inMi+3 = {EJ rG[o,T]) namely, the Schwarzschild orthonormal frame (|3.10l); 
Figurej^ Let {eijg, eg = do, be a frame extension in expressed in terms of 

dd via 


(4.1) 


= O^dd 


Consider now the (unique) metric g for which is orthonormal, gab ■= "m-ab = 
diag(— 1 , 1 , 1 , 1 ), and the corresponding connection coefficients (Ay)^^ = g{Ve^ei, ej). 
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Then Proposition 2.1 asserts that the EVE for g, under the Lorentz gauge 

conditioiH 

(4.2) (div^x) = (^3)23(^2)ij, 
reduce to the system of scalar equations 

ek{Af^)ij + A^])ij + e^((A[^)^] (Ak)ij) 

(4.3) + + {[A>^,Ffj_^])ij + e^({A3)23iA2)ij) + e„({A^)fj.''{Ak)ij) 

^o(Oc) = ~ {A[Q)h]‘^ + (A[o)c]^0^, ly, i,j, c,d G {0,1,2,3} 

where □ := -Cq + + e§ and '®(A[o)b]‘^ = [5o,9t,]‘^. 


However, the system (4.3) has singular initial data in the Schwarzschild background 
which do not permit an energy approach directly. Eor this reason we recast the 
equations in a way that captures the closeness to the Schwarzschild spacetime. Let 

(4.4) (iijy)ij := {A^)ij — ^(Ai,)ij : {ET-}.rG[o.T] “t M GJ G {0,1,2,3}, 

where the components ^(A^)i -i are the Schwarzschild conn ection coefficients corre¬ 
sponding to the frame {9i}o (3.10) and they are given by (3.11). We are going to 
use these new functions to control the evolution of the perturbed spacetime. 


Consider now the analogous system to (4.3) sa tisfied by the Schwarzschild compo¬ 
nents ^{A,j)ij,dc. In view of the asymptotics (3.15), we define P, to be a smooth 
function satisfying the bound 

Cq,k 


(4.5) 






j.q+^k 


for constants Cq^Cq^k depending on M >0. Taking the difference of the two 
analogous systems we obtain a new system for the functions — written 

schematically in the form: 

h’^^dadbMij =OT:idu + OT^u -f 0r|(0 - /) -f OT^diO - I) 

(4.6) +T 3 V? + Oudu + 0d{0 — I)du 

do{Oi - Ic"^) = Pa (O - /) + (O - I)u + u, 

where 

(4.7) := = 5 “'’ 

and each term in the RHS denotes some algebraic combination of finite number of 
terms of the depicted type (varying in J,j) where the particular indices do not 
matter. 


Remark 4.1. Evidently, the systems (4.3) and (4.6) are equivalent. The benefit 
is that the assumption on the perturbed spacetime, being close to Schwarzschild, 
implies that the functions — are now small and regular. Thus, we have 

reduced the evolutionary problem to solving the PDE-ODE system of equations 
(4.6). However, the issue of singular initial data in (4.3) has become an issue of 
singularities in the coefficients of the resulting equations (4.6), at r = a: = 0, which 
do not make it possible to apply the energy procedure in standard spaces; see also 


12- 


We choose now a specific type based on the one satisfied by the Schwarzschild reference 


frame i 3.14 i. 
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(3.17). These singularities, in large part, are due to the intrinsic curvature blow up 


and cannot be gauged away; in particular the coefficients r 3 of the potential terms 


in (4.6) correspond to the Schwarzschild curvature (1.4). Some of the functions T 


that appear in (4.6), expressed in terms of Schwarzschild connection coefficients 


(3.11) and their derivatives, are less singular than (4.5), but representatives of the 


exact bound do appear in all the terms. 

Remark 4.2. Another crucial asymptotic behaviour that our method heavily de¬ 


pends on is that of the radius function r. According to (3.6), we observe that the 
best bound one could hope for the ratio l/r^ is of the form 


(4.8) 


, 1 „ C 


which obviously fails to be integrable in time r S [0,T], for any T > 0. This fact 
lies at the heart of the difficulty of closing a Gronwall type estimate. 


4.2. The weighted spaces. In order to study the well-posedness of (4.6) we 


introduce certain weighted norms. It turns out that the weights which yield the 
desired flexibility in obtaining energy estimates are the following. 

Definition 4.3. Given a > 0 and r G [0,T], we define the (time dependent) 
weighted Sobolev space as a subspace of the standard iJ® space on 

with the Schwarzschild induced volume form satisfying: 


(4.9) R®’“[r] : w S iJ®(S^), 




-E 


k<s 


'Sr 




where by 9*^^^ we denote any order k combination of directional derivatives with 


respect to the components 81 , 82,83 of the Schwarzschild frame (3.10). For conve¬ 
nience, we drop r from the notation whenever the context is clear. 


Remark 4.4. Observe that the weights in the norm || • in (4.9) blow up only 


at T = 0, a; = 0. For r > 0 fixed, the weights are uniformly bounded above by some 
positive constant Cr, which becomes infinite as r —)■ 0"^. The dependence of the 
power 2a — 3(fc — 1) on the number k of derivatives corresponds to the singularities 


in the coefficients of the equation (4.6). 


Lemma 4.5. The weighted spaces satisfy the properties: 

R®i’“ C Si < S2 

(4.10) r~i’-u G whenever u G i7®’“ 

aWuG 77®-'=’“-?'= k<s,uGH‘^’‘^ 


Proof. They are immediate consequences of Definition |4.3| and and the fact that 
|9i(r-i')| < Clr-i^--^ \d2{r-"^^)\ = |93(r-5')| = 0, 


cf. (3.5), (3.10). 


□ 
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4.3. Local existence theorems. Let 

3 

sup (||K)*j||ff3,„ + \\do{u^)ijf 

.rG[0,r] ^ 


£{u,0-a,T):= ^ 

pT 


.-4 

1 o / 


(4.11) 


{WMi. 




3 r 'X' 

^ sup + [ wot-I, 

LrGfO.Tl ^ = Jo 


d ||2 

11^3 


i dr 


c , d =0 L ^6 [ 0.^1 

be the total weighted energy of the functions {u^)ij,Ot — Ic^ defined in {Ei-},-g[o.T] 
(3.1), Figure 3.1 the backward domain of dependence of Eq with respect to the 
metric g we are solving for. Since the actual domain depends on the unknown 
solution, it will be fully determined in the end; see Section For brevity we 
denote by 


(4.12) 


£o ■— 

i>',2,jG{0,1)2,3} 

+ E wot-I 


u^)ijiT = 0)||^3.c + WdoMiji^ = 0)11 




d \\2 

11^3 


c.dG{0,l,2.3} 

the energy at the initial singular slice Eq. 


■2(i:o) 


Remark 4.6. It would seem more natural if the components {uu)ij were lying 
in that is, one derivative less than Ot — Ic'^- However, this would create 

additional technical difhculties, which we choose to avoid by treating both 
ot — Ic‘^ at the same footing (number of derivatives), since it is an issue of the 
structure of the system (4.6) and not a singularity issue. 


The following theorem is our first main local well-posedness result for the system 
(4.6), whose proof occupies Section 

Theorem 4.7. There exist a > 0 sufficiently large and £ > 0 small such that if 

(4.13) £o<+oo <£, c,d = 0,l,2,3, 

then the system (f.d) admits a unique solution, up to some small time T = T(£q, a) > 
0, in the spaces 

G C([0,r];R3.a)nL2Qo,T];R3.a+i) ^^0 £ {0,1, 2, 3} 

(4.14) ao(u.)., £ C([0,r];R2.«-f)nL2Qo,T];R2.«-{) 


Ot - 1/ £ C([0, T]; ) p L^{[0^ T]; R^.a+I) 


c,de 10,1,2,3} 


Remark 4.8. (i) The second part of condition ( 4.13| ), £ > 0 small, is necessary 
for th e eq uation (4.6) to be hyperbolic, yielding sufficient pointwise control on the 
/i“'”s ^ 

1 
2 


(4.15) 




< Ce^ 


6, c = 0,1,2,3, bfc. 


It could be obviously replaced by the stronger assumption that Sq < e, since the 
energy £{u, O; a, T) controls the L°° norm of u, O by standard Sobolev embedding. 
{ii) How large the exponent a has to be depends on the coefficients of the system 
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(4.6). In the final inequalities in ^ a > 0 is picked large enough so that certain 
‘critical’ terms can be absorbed in the LHS and the estimates can close. 


The above theorem is a local well-posedness result for the system (4.6). Imposing 


now the proper conditions on the initial data set of (4.6), the solution (4.14) yields a 
solution of (4.3) which in turn corresponds to an Einstein vacuum spacetime (O- 


Theorem 4.9. Let a,e be such as in Theorem \ 4.1 and let (T,Q,g,K) be an initial 
data set for the Einstein vacuum equations {If satisfying the constraints 
such that the components 


(4.16) 


v,i,j = 1,2,3, 


(4.17) 0,^-/e‘'ei73.«H-|(Eo) ||0f-4‘'|U=o(s„) <e c,d=l,2,3, 

computed with respect to an orthonormal frame {ei}J on (Eo,g), and 

(4.18) (n,)oj(T = 0) Gi73-“(Eo) z,j = 1,2,3. 


Then, there exists a solution g to the EVE (1.1) in the backward region to Eq, foli¬ 


ated by {ET}tg[o,T]; with induced initial data set {^,K) on Eq and an orthonormal 
fra me e xt ensi on {e^jQ for which the co rresponding (spacetime) functions {u,f)ij, Of:— 
if (^.j)),H.l) lie in the spaces (f-l)^- 


If in addition Of — Ic'^ G (^([O, T]; )^ c,d = 1,2,3, then the Einsteinian 

vacuum development is unique up to isometry. 

The fact that such (non-spherically symmetric) initial data sets (T,Q,g,K) exist, in 
compliance with Theorem |4.9| is shown in 

Proof of Theorem\4.9\ We want to invoke Theorem |4.7[ For this purpose, we pre¬ 


scribe initial data for the system (4.6): 


(i) The components (4.16), ( 4.17 ), (4.18) are given. 

(ii) Since in the beginning of (4.1 we assumed cq = do and since is initially 
tangent to Eq, we set 

Of{T = 0)=Ia° 


(4.19) Of{T = 0)=Io 
(in) We (freely) assigij^ 

(4.20) (uo)ab(r = 0) := iAo)ab - HAo)ab G i73’“(Eo), 


a, 6 = 0,1, 2, 3. 


a,b = 0, 1, 2, 3. 


Once we have prescribed the above, the components do(uv)ij(T = 0) are fixed 
by the assumption (2.21) on the initial data of the original system ( |4.3[ ); see Re¬ 
mark |Ai2 Indeed, subtracting the corresponding Schwarzschild components from 


(A.15),(A.16), which obvisouly satisfy the same initial relations, cf. (3.14), we 


obtain schematically: 

(4.21) do{u,),j = Odau + Tiu + r3(0 - /) + ' 


on Eq, 0 = 1 , 2 ,3. 


By (|4.10|) and standard Sobolev embedding we conclude that 

^ Obj = 0,1,2,3. 


(4.22) 


doMijir = 0) G 


^^The functions (uo) „h(T = 0) or equiv alent ly (Aolatlr = 0) fix the do derivative of the frame 
{eilp on So; see Lemma A.l and Remark A.2 
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Thus, the assumption (4.131 is verified and Theorem 4.7 can be invoked. From 
Proposition 2.1 it follows that the solution ( |4.14| ) of (4.6 1 and hence of (4.3) yields 
indeed an Einstein vacuum spacetime ({St- }G[o.r],5)- 


To pr ove uniqueness (up to isometry) we rely on the uniqueness statement in The¬ 
orem 4.7 Suppose there is another Einsteinian vacuum development of 

the initial data set {Tto,g,K), diffeomorphic to {ET-}T-g[o_T]; satisfying the hypothe¬ 


sis (4.16), (4.17), (4.18); defined by pulling back the relevant quantities through the 
preceding diffeomorphism, taking differences etc. In order to use the uniqueness 
statement in Theorem \A.7\ we need the two spacetimes to have the same initial 


data for the system (4.6). The part of the initial data set given by the assumptions 

The 


in the statement of Theorem |4.9| is of course identical for both spacetimes, 
remaining components that we want to agree, other than the (uo)ab(T = 0)’s, as 
noted in the previous paragraph, can be fixed by condition (2.21). Therefore, we 


get identical initial data components for the system (4.6) by constructing a Lorentz 
gauge frame (4.2) {ei}g for g, which is initially equal to (eilg on Sq and such that 


({(o)ob('F = 0) = {uQ)ab(T = 0) as well; see Lemma A.l The only assumption to be 


verifie d is the well-posedness of the system (A.l) for functions in the solution spaces 
(4.14), after taking differences with the equation for the frame {e^lg. However, this 


falls in the category of the system (4.6) [in fact simpler, being semilinear] to which 
Theorem |4.7| can be applied. The extra derivative that we have to as sume in order 
to close, of — Ic‘^ S is due to the divA term in the RHS of (A.l). □ 


5. Proof of Theorem 14.71 

Throughout this section we will use the notation X < Y to denote an inequality 
between the quantities X, Y of the form X < CY, where C is an absolute pos¬ 
itive constant depending only on the Schwarzschild mass M > 0. The same for 
the standard notation 0{X), for a quantity bounded by |0(A')| < CX, X > 0. 
Furthermore, all the estimates regard only the Schwarzschild region foliated by 
{Sr}rG[o,T]; Figure]^ 


5.1. Proof outline. We prove Theorem |4.7| via a contraction mapping argument. 


First we establish an energy estimate in the relevant weighted spaces in [5.3 


Then we obtain a contraction, in [ 5.4 in the corresponding spaces of one derivative 
less, see (5.47), which together with the energy estimate yield the desired solution 

( |4T^ . 

To derive these estimates we have to eliminate some critical terms which are gen¬ 
erated due to the singularities in the coefficients of the equations, having larger 
weights than the ones in the norm (4.9), and which prevent us from closing (see 
Propositions 5.2|5.5 ). This is where the role of the weights (4.9) comes in. The 
parameter a > 0 helps generate critical terms with a favourable sign. Being large 
enough, but finite, a provides an overall negative sign for the critical terms, hence, 
rendering them removable from the RHS of the final inequalities. This enables us 
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to close the estimates and complete the proof. The precise asymptotics of the sin¬ 
gularities in the coefficients of the equations (4.6 1 , at r = cc = 0, and the opening 
up rate of the radius function r in r > 0 play a crucial role herej^ 


5.2. Basic estimates. Let u be a scalar function defined on represented by 

(5.1) V o ipT- : Ur ^ 

where ipr ■ Ur ^ E,- is the {x,d,(p) coordinate chart. We recall some standard 
inequalities: the classical Sobolev embedding of in L°°{U) and the interpo¬ 

lation inequality 

(5.2) MLHu)<C\\v\\l^^^\\yv\\l^^^, 

for a bounded domain U C with (piecewise) boundary. In the following 
proposition v is assumed to be regular enough such that the RHSs make sense. 

Proposition 5.1. For a general function u : E,- — )■ K, r G [0,T], with the appro¬ 
priate regularity, the following inequalities hold: 

The L°° bound 

(5.3) + 

and the estimate 

(5.4) + 


Proof. From the embedding H‘^{Ur) ^ L°°{Ur) we have 


(substituting (3.8) and the frame (3.10)) 


< (k+iy 






We argue similarly in the case of (5.4). 


□ 


5.3. Energy estimate in We define now the mapping, which really corre¬ 
sponds to an iterative process. Let {u, 0 | := |(M^)i,-, : v, i,j, c,d = 0, 1, 2,3} be 

a set of functions in the solution spaces (4.14), verifying \0^ — Ic^l < e initially on 
Eq. Let T > 0 be sufficiently small such that 


(5.5) 

We also assum^Hl 




(5.6) ||ao(0:)f 


£{u^ O; a, T) < 2£q. 

£q -\- £o Vt € [0, T], c, d = 0,1, 2, 3. 


we were to tweak the leading orders just by e > 0, the previous procedure would fail no 
matter how large o > 0 is to begin with. 

^^Any assumptions that we make on the functions u, O, we must derive for the next set of 
functions u, / below. 
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Iteration step'. Consider the following linear version of the system (4.6), where we 
replace the functions m, O in the following specific terms by the corresponding ones 
from the set {u, O}: 


(5.7) 


h'^^dadb{u^)ij = OT sdu + OT^u + OF 9(0 — /) + 0 + 30(0 — I) 
+ T^TL^ + Oudu + + 00(0 — I)du 

do{Oi - = Fa (O - /) + (O - I)u + u, 


-^ab 


-^a—=rb 


where h = Observe that we kept in the RFIS of (5.7) the functions 

u, O attached to the most singular coefficients of the system. This is actually very 
important to our strategy. Flad we replaced them with the corresponding functions 
M, O as well, it would not be plausible to derive a weighted energy estimate. 


We assume now there exists a solution (uy)y, 0(f — Ic of (5.7) lying in the solution 
space (4.14). The existence of such a solution is based mainly on the energy estimate 
we will derive below and a standard duality argument which we omit. 

Claim'. For a chosen large enough a > 0 and T > 0 sufficiently small (depending 
on Eq, a) the following estimate holds 


(5.8) 


E{u, O; a, T) < 2£o- 


The preceding //^-weighted energy estimate, cf. (4.11), will be used in the next sub¬ 
section to close the contraction argument that yields the existence and uniqueness 
of the solution (4.14) to (4.6). Now we begin the proof of (5.8): 

First note that by the fundamental theorem of calculus, following a 0o integral curve 
and employing (5.3), we readily obtain from our initial assumptions and (5.6) the 
pointwise bound 

(5.9) 


sup 

tG[0,T] 


\0 — /|||l“(s^) < £ + CTEq < 2e, 


provided a > ^ + 3 + ^ and T < . 


All the more , direc tly from the ODE in (5.7) we deduce the esti mate: [applying the 
bounds (5.3), (5.5) to [O — I)u and employing the asymptotics (4.5)] 


( 5 . 10 ) + 

for all re [0,T], c, d = 0,1, 2, 3{3 


Ml.. 


[r]> 


We derive (|5.8|) in the backward domain of dependence of Eq w.r.t. the metric 


-^ab 


{gab)u ■= ff«(0a, 0&), a, 5 = 0,1, 2, 3, whose inverse is given by := h ; compare to 
(4.7). The boundary of the domain is the backward incoming gn-null hypersurface 
emanating from 0Eo (Figure]^. We foliate the domain by the r = const, 
hypersurfaces inside . Let p be the scalar function defined near Af” via 


(5.11) 


p(C“) :=r-r, 


16. 


This estimate, together with (5.8i in the end, imply the analogue of (5.6( for the functions 


doioi). 
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where C“ is the cylinder obtained from the flow of 9S“ backwards along the integral 
curves of Oq. Using p we may write each leaf of the foliation as 

(5.12) U {p,=T-t*}\jBr tG[0,T], 

fG[T,T] 

where Pr '■= p\^u is simply the projection of onto S“ through the 

integral curves of dg. Since by definition p + t — T is zero on it follows that 



Figure 6. 


the (iij-gradient of p + r — T, on lies on the hypersurface itself and furthermore 
it is p^i-null, i.e., p satisfies the eikonal equation 

(5.13) \Vg^{p + T - r)|g^ = dA{p)dB{p) + + 211-1/1^°5a( p) 

= 0 on 

where A,B = \, 2,3. 


In this regard, we define the following adapted energy, which controls the part of 
the total energy (4.11) that refers to u: 


(5.14) 


Es+i,a[u]{T) :=i Y 


^00 [5o(Wi/)ij^j] 


j,2a—3\J\ 


^ ^AB dA{p^u')ij^ 

' R ^ 37Ti I sm “T 


{n.? 


ijj 


^-lui 


.-lui 


^2a-3(|J|-l) 


dps 


where {u^)ij^j := 5i'^i(u^)ij and J is a spatial multi-index (containing only direc- 


-^ab 


tions 9i, 92 , 93 ). It is evident from (5.9), h = that is equivalent 

to the weighted x '2 norm of u on E“. 


We summarize in the following proposition the main energy estimates derived below. 
Proposition 5.2. The following two energy estimates hold: 

drEs^alu] + 8M^e-^(l - r)a£'3,„+iM 
(5-15) ^ i^O + ^0 + + Q^^'^'o)^'3,aM + Fls.Q-l-lM + fo||0 —/|P 3 ^^3 

+ \\0 - 4||^3 ,^^5 + a^£o + Sq 
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(5.16) 


\\Oi - + 4M2e-i(l - T)aY,\\Oi 

c^d c,d 

^ \\0 - 4||^3 + £'3,a+lM + £l, 


d ||2 

c 


,a + 


5 

2 


for all T G (0, T). 


The overall energy estimate (5.8) follows from Proposition 5.2 Adding (5.15), 


(5.16) we wish to close the estimate by employing the standard Gronwall lemma. 
However, this is not possible in general, because of the critical energies in the RHS, 
having larger weights than the ones differentiated in the LHS, namely, 

\\0 — /|p 3 ^ 5 instead of i? 3 _Q,[u], \\0 — I\\^ ^ ^ 3 . It is precisely at this point that 
the role of the weights we introduced is revealed. Choosing a > 0 large enough 
to begin with, how large depending on the constants in the above inequalities, we 
absorb the critical terms 


A3,„+iM,||o-j|!^3, 


a + 


5 

2 


in the LHS and then the standard Gronwall lemma applies to give (5.8). 


Proof of { 5.15). Let 


(5.17) := - 


— h 


00 


[9o('U^) 




j.2a—3\J\ 


H" ^ ^ irm ^ irm r 




ij.J 






rp2cx-\-3 —3|^| 


for any spatial multi-index J with |J| < 2; recall {uv)ij,j '■= d^^\ui,)ij. It follows 
from (5.12) and the coarea formula that 


(5.18) 



/as; 


Pj,a 


dS + 


dlrPj^ad^Sg 


P.J.a^rdfiS g ^ 


where stands for the gradient of p with respect to the intrinsic connection on 
{Y,ri^g) and dS is the S chwarzschild induced volume form on 9E“. Note that the 
boundary term in (5.18) has a favourable sign. Since A/^ is ^ij-incoming null, we 
show that the sum of all arising boundary terms has a good sign and therefore can 
be dropped in the end. To analyse the last two terms in (5.18), we recall the d^. 
differentiation formulas of the radius function r (3.5), the estimate on volume form 
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d^Sg (3.9) and the commutator relation [9 o,i9b] = ‘^(^[o)b]'^5c 

/ d)rPj^adfJ,Sq I Pj^a^rd^Sg 

JS" 

= -8M^(l-T)a/ e~™Pj^a+idfJ.sg+ [ Pj^aO{^)diisg 
JS" JE“ 

T-OO. „ fj-AB dA{u^)ij,J dB{u^)ij,J 


(5.19) 






- +ao(/. ) 

j^O dQ{Ui/^ij^jdQ(Ui/^ij^J ^ -j^AB dA{j^u)ij^J (^1^)2^, J 


^2a —3| J| 


^cn—||J| ^a—||J| 


djJLSg 

d/iSg 


J-AB dA{Uu)ij,J 


f^iUv)ij,jdo{Uu)ij,J J 


The first term on the LHS of (5.19) is critical having a favourable sign of magnitude 
a. We use this term alon e to absorb all arising c ritic al terms in the pro cess . Recall 
1^1 = |0 I < 1, cf. (5.9), and the asymptotics (4.5). Also, applying (5.3) to Ogh 
and (|5.6l) we derive 


IQdo(h)l<£(u,0;a,T)2, 


1 


7* 2 


|r3|<4- 

7* 2 


Hence, by Cauchy’s inequality and ( |5.5[ ) we have 


+ / £lh 


c, ij:°^A9o{Uu)ij,j\ ^ fj-AB dA{u^)ij,J dB{u„)ijJ 
) ^2a-3|J| +^0W ) 

-rABdA{u^)ij,J r|i9(Uj.)y,j 


||J| ||J| 


d^Sg 


2 


-d^s 


_ ^ ^2a+3-|J| 


J 7* 2l I 7* 2 I 

(5.20) < £■ (m, O; a, T) 5 [u] + [u] 

;< ■S’d'^S.aH + A3^a+i[M] 

For the next term we proceed by integrating by partj^ (IBP), denoting by A^ := 
NbOb the outward unit normal on 9E“ w.r.t. Schwarzschild metric g on S“: 


n 


j^O Ogi^Ui^^ij^jdgiUi^^ij^j ^ j^AB dA{^Ui/)ij,J dBdgi^Uii'jij^j 


T"*-*- 2 I I 7* 2 I 


d^s 






j^oo d‘^{u^),jj -^AB dBdA{u^)ij,J 


j,a 2 I "tl 


a-f |J| 


d/is 


/as; 


- 2 I >^1 2 I 




^s(j,2q-3 |J| )^A{Uu)ij,J + 


—r 3 j 


l^2a —3| J| 


^0 {^^l/)ij,jdlJ-S g 


integrate by parts using the spatial part of the Schwarzschild frame ^1,52,53- Doing so 
we pick up connection coefficients, since it is not covariant IBP. 
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It is immediate from the definition of the frame (3.10) and (3.5) that 


In ^ ^ M ^ “ 

2fv-.s) ^ „2fv- 


^2(^) = ds(^) = 0. 


2 


7*^ 2 


Hence, similarly to (5.20) 


(5.22) 


< j.2a — 3\ J\ 


) J “1“ 




j.2a—3\J\ 


^0 )2_7, g 


(I J| < 2) < (£| + a^)E3Ay] + E3,a+i[u]. 


Remark: The term in the RHS of the preceding estimate with coefficient is not 
critical. This is very important otherwise the overall estimates would not close, 
since the critical term with favourable sign in (5.19) is only of magnitude a. 


We proc eed to the boundary term in the RHS of (5.21). Recall that p is constant 
on i9E“ (5.11), and decreasing in the interior direction of E“. Hence, the outward 
unit normal N is the Schwarzschild normalized gradient of p on E“, N = . 


Since ih )a,b=i,2.3 is a symmetric positive definite matrix, the following standard 
inequality holds: 


j-ABdA{Uy)ij^j 


< 


(by (5.13)) 




(n^h^'^NANB) 

_4 R 

n^h dA{p)dB{p) 

|‘5Vp|2 

-h°° - 2nh^°dA{p) 

|^Vp|2 


Therefore, we have the bound 


(5.23) 


^j-AB dA^jy j 


IdTA 


.-iui 


.-iui 


< 



^0 


ij,J I 


.-fUl 




-h°° - 2 nh^°dA{p) 


I^Vpl 


^AB 


\ I^Vpl ■ 


^A{^ 


)ij,J 


dBiUi,) 


I' )ij,J 




.-lui 


dS 


< 1 /■ [go(».).j,j]" 2nh^°dAip)[doM,j,jf 

- 2 igs- I^Vp^l r2“-3|J| \^Vpr\ r2«-3|J| 

_ 
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The remaining term to be estimated is the one on first line in the RHS of (5.21), 
which we rewrite 


n 




T 2 


-ii^i 


7*^ 2 


-flJI 


^0 


T 2 


-fl^l 


d^s 


(5.24) =-j 

JE“ ^ 

f j-AodAdo{u„)^jjdo{u^)^jj r3i9('u^)y,jao('u^)*j,j 

+ 2ff/l -^ + nh^ -- dfiSg 

By taking the derivative (J spatial multi-index | J| < 2) of the first equation 
in (5.7) and commuting the differentiation in the LHS we obtain the equation 

Tl'^dadb{u^)ij,J 


(5.25) 




OTsdu + OTsu + OTi {O - I) + OTsdiO - I) 


+ + Oudu -I- -I- 0d{0 — I)dv 


[h‘'^dadb,d^'^'>]{u„)ij, 


where the commutator can in turn be written schematically as: [recall (3.15 ),(4.5)] 


(5.26) [h°'^dadb,d^'^^]{u^)ij = d‘^{h)d'^{u^),j + [T3d{h)+ r3h]d'^{u^),j 

+ d(h)d^(u^)ij + hTid^{u„)ij 
+ [d(h)T3 + hT9_]d{u,),, if|J|=2 

[li'^dadb,d^'^'>]iu^)ij = d(h)d‘^{u„)ij +hT3^‘^{u^)^J if jJj = 1 

+ hT3d{u^)ij 


We integrate by parts in the second term on the RHS of (5.24) and argue similarly 
to (5.22) to get 


2nh 


(5.27) = 


E“ 


'3e: 


■p 2 cic— 3 \J\ 


^2l0 [doiUi^)ij,j] f (•Tr^0\ [^o{Ui^)ij,j] J 

ndAih ) ■ 

( ^2^%J\ ) [^0 j] ^dflSg + fnhT I dfiSg 

—r3 d{'Ui/')ij jdoi^iii/^ij J 


+ / rth 


< 


flh 


JdY.1 


+ C(4 + a^)E^An] + 
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Finally, for the last and main term in the first line of the RHS of (5.24) we recall 
that |n| < ^ to obtain directly from Cauchy’s inequality 

rl 

- [_ {h‘'^dadb{u^),jj)n dfisy, 


(5.28) 


< 


< 


I h"'’dadb{u^)ijj ||2 
dadb{u^)ijj u2 


j,2a—Z\J\ 

I ^a—I |J| + 1 "L^ 


li ^0 




We proceed by plugging the RHS of (5.25) into the first term in the last inequality 
(5.28) above and treat eac h ari sing group of terms separately. Empl oyin g the basic 


inequalities in Proposition 5.1 along with the bounds of 0,d{0),u (5.5), (5.6) and 
(5.9) we derive: 


(5.29) 


[Oraan + OPaM] 


lL2 


21 "^1 2 

< |io|!i..F;3.c.+iM + \\dm\l<^E,,M + \MU. 
d(E(p)T.du d(E(p)T,u , 

+ II--IIl2 + II-^^.iiJiji-'IIl^ 

(the last two terms appear only in the case | J| =2) 


< £{ u , 0 ] a , T ) E 3 ^ a [ u \ + Es ^ a+lM + (||m 
^ ^0^'3,q:M + -E'a.a+lM 


+ Il5«lli==)ll^^lli^ 


(5.30) 


[ora (o - /) + or3a(o - /)] 


pOt 2(^1 2 


lL2 


<110-^11 


+ 


H->,- 

, 0 - 1 . 


f 2 


.^j + liao||i.||o-/|| 
i» + lia(0-j)||i„ 






(we include the last term only when | J| = 2 and utilize (5.6),(5.10)) 


<l|0-/f ,o+s+(^o+^o)||0-- 


H' 




+ £o||0-/|| 


H-’ 


(5.31) 


+ + £o||M||^2.a 

I dO'> [Fa?!^ + Oudu + # + 0d{0 — I)dv\ ,,2 


< 


j,a. ^\J\ 2 

l4lli- + ll?lli- + liao|ii. 

7* 2 r 2 


lL2 


hilio 


liauiii. 


E3,a\u\ 


Nii„||ah||io„||o-/|| 




2,„+j -F tll“llL' 

d^O-I)d^u 


|iz||i»l|ah||i»)ll 


M\h.. 


„a—3— k 


\\do\\UE3,4u] + 

(liao||i» + liao||i=„liah||io„)(l|o-j|| 


i2 + ||aii||i„||o-/|| 


H-’ 




_3 + E2^a[u] 
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^2 , , II '5^(0) ||2 II ||2 


< E{u, O; a, Ty + £{% O; a, TY + 


Tlli^ll- 


T 2 2 4 


rlli^ 


(employing the estimate (5.4)) 

<Sl+£l + aHl 

By (5.29)-( 5.31[ ) we have the following lemma. 

Lemma 5.3. dQ{uv)ij G (^([O, T]; n L^([0, T]; and moreover the 

following estimate holds: 

II IL=[t] 

(5.32) < £’o(£'3,aM + \\0 - /||^3_„+3 ) + Es^a+lM + ||0 - 


/or|J|<l, JC{1,2,3}, rG(0,T). 


Proof. The proof follows by solving for d^iu^ij in the equation (5.7) and summing 
up the above estimates (5.29l-(5.311. □ 


To bound the commutator (5.26) we treat the cases | J| = 2, | J| = 1 separately. For 
|J| = 1: 


(5.33) 




Ck_3_ 1 
7 * 2 2 


IL2 


d{h)d'^{ui,)ij + hT3d'^{ui,)tj +hT3d{u, 


a_3_ 1 
T 2 2 


||2 
Il2 


< 


liahii 


2 II ^ II 2 


^Q: r. o 


^Ili3 


I 2 II ^ {uYij II 2 


Flli3 + 


lio 


I II 2 


„a+l 


Ili3 


(employing (5.32) in the case d'^luYij = '9oK)*j) 


^ ^o{E 3,a[u\ + ||0 - /||^3_„+3 ) + if3,a + lM + \\0 - /|i^3 
+ a3£:2 + £3 


-0 ^0 
When I J| = 2 we have 


[/l“^l9al9b,a(-^)]('U^)y ||2 


(5.34) 


y.OL — 3 — 7 T 


lL2 


< 


d^{h)d^{ui,)ij + [r3 5(/i) + r3h]d‘^{ui,] 


^a— 3 —7 


U l|2 


lL2 


I d{h)d^{uyij + hVid^iuYij + [9(/i)r3 + /ir|]5(n^)jj 2 

IIl2 


— 3— o 


(note that term d^(ui,)ij contains at most two do derivatives) 


92/1 


2 II ||2 I ||q7:i|2 / 77 L,1 I 1102/ ' "2 


13 ||//l,<7-3; 


< Il^^^lli4||^^^lli4 + \\dh\\iyE3An] + WdlAY 

7*224 7*224 

+ ll^lli“ {E3,a+l[u] + ||9o(Uj7)y 11 ^ 1 , 0 - 2 ) 

(employing the A estimate (5.4) and (5.32)) 

< AEqEoAu] + £oE3A^\ + -E'S.a+iM + folio — 7|P 3 „^3 + ||0 — 7|P 3 

rz ’ 2 7i 
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cx £q £q 


Summary: Incorporating (5.19)-(5.34) in (5.181 we conclude that 


(5.35) 


dr Pj^adfJ-Sg+ 8M^e ^(l-T)a/ Pj,a+idnsg 

^ (^0 + + i?3,Q+i[M] + ||0 — /|| 




+ £o\\0-I\\ 




+ a £q +£q 


Summing over the indices i,j and J, |J| < 2, we arrive at the desired estimate 

([51^. □ 


Proof of {5.16). Let J, | J| < 3, be a spatial multi-index. Like in the case of (5.18), 
it follows from the coarea formula and the asymptotics (3.5),(3.9) that 


(5.36) -a 


1 „ „ Of, - I 


,J ~ ,J ||2 




2“'" ^a+3-§|J| 2 |S'Vp|r2“+6-3|J| 


dS 


3|JI f 


- (a -I- 3- —) 


j»2q;+ 7—3|t7| 


■drrd^s 


+ ^2a+6-3lj| 




® ' 2 j.2a-|:6-3iJ| -drdfiSg 


< -AM\ 


+ / n 


{ojj - i^\j)Moij) 


dlJ,Sg + C^ll " 


— T 


,J ^ ||2 


j.2a+6—3\J\ I ^11 ^^_|_4_3|j| NL2(S^)’ 

where O^ j - := - l/). By Cauchy’s inequality we have (^2 y) 


(5.37) 


n 


iot - h\j)d^{oij) 


■Oc,J - Ic ,J 2 


- 1 . ,U /-yjK- C,J / 1 ^ II “ c,a ||2 

r2a+6-3|J| - II ^a+4-S|J| 


Moij) „2 

„+s_S|j| Wl^ 


Taking the derivative of the ODE in (5.7) we obtain 
(5.38) doiO^j - 1/j) = d^-^'> r3{0-I) + (0-I)u + u 


+ [do,d^^^]{Oi-l/) 


The commutator in the RHS of (5.38) schematically reads 


= T3diOf-l/) 

if 1 Jl = 1 

= TsdiOi - I/) + T3d^{0^, - I/) 

if 1 J| = 2 

= T 9 d{oi -1/) + r^d^oi -1/) 

if 1 J| = 3 

+ T3d^iOt-l/), 
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where we note that at most one do derivative of — Ic‘^ appears in the preceding 
expres; 

(5.40) 


expressions. Hence, we deduce directly from (5.38): 


'Ih 


^ l|r|(0 — /)||^3,c» + l + 11(0 — /)u||^3,<. + l + +11^11^3,<» + l 
+ 11 ^C+I-IUI 


(exmploying Lemma 4.5 and applying the bound on (O — I)u) 


<llo-j|it3,„+a+ l;3,.+iM + I|o-/|P 


'H 


3,a+^ 

'\d 


'H 

d\ 


3,a + 


3||^||//3,< 


T^d[oi - //) r39(o,^ - 4'^) + T^d^ioi -1/) 

' '\h + \\ -- Wh 


^a+|—3 

T9d{of -1/) + -1/) + r3d^{o^-1/) 


j,Lj.-r 2 2 T 

-\d T d\ , T_;:)2//nd t d\ , -n„^i3//0(i T d\ 


r^~ 2 2 


2 

L2 


^ l|0-^f^3,.+s + £{%0-a,Tf + Eo,o.+i[u] 


Combining ( |5.36[ )-( |5.40[ ) we derive 
(5.41) 




f^d T d 

I ^c,J - -'c ,J ||2 


j.a-\-Z— 2 01 

<||0-/f 


H 


3.Q + 


i2(£u)+4M e q;( 1 —r) 
5 + E?,,a+i[u] + £o 


f^d T d 

.J||2 

a+4-||J| IIlHS?) 


Taking into account the set of indices c, d and J, | J| < 3, we complete the proof of 
(5.16) and hence of Proposition 5.2 □ 


5.4. Contraction mapping in We proceed to show that the mapping de¬ 

fined via (5.7) in the beginning of (5.3 is a contraction. Let us consider another 
solution set {u^)ij,0^,u,0 of the analogous coupled system (5.7). Setting 

(5.42) {du^)ij = {uv)ii — {uv)ij, du = u — u, dO'!: = Ot — Ot, dO = O — O 


we obtain schematically the new system of equations (depicting only the types of 
terms in the RHS suppressing the particular indices) 

h°'^dadb{du^)ij 

(5.43) =OT:id{du)+ Or3du + OT3dO 

+ dO[T3du + Tou -k Pa(d - /) -k Ps^O] 

-k OTsdidf) + (0 + d)dOd^{u^)ij + G{du, dO), 

where 

G{du, dO) = Ta du{u -I- u -I- -I- -I- uu) + Oud{du) 

(5.44) -I- Odudu + dOudu + 0d{0)d{du) 

+ Od{dO)du + d0d{0)du 
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and 

(5.45) 9o((iO)?) = r adO + (O — I)du + udO + du 

Further, we assume that both sets of variables we have introduced are consistent 


with the energy estimate (5.8) we have established in the previous subsection: 

(5.46) £{u,0;a,T), £{u,0]a,T), £{u,0\a,T), £{u,0\a,T) <2 £q. 

Claim: For large enough a > 0 and T > 0 is sulRciently small the following 
contraction holds: 


(5.47) 


'H 


2,a+| 7 5 


E2Adu] + < n{E2Adu] + \\dO\\] 

c,d 

for some 0 < k < 1. 

Remark 5.4. We are forced to close the contraction mapping argument in 


having one derivative less than the space of the energy estimate (5.461, see (5.3 as it 
is common in 2"'^-order quasilinear hyperbolic PDF because of the problematic 


term (O + 0)dOd^{uu)ij in (5.43), which is generated from the difference of the 
top order terms in the LHS. 

Proposition 5.5. Under the above considerations, the following estimates hold: 
(5.48) drE2Ad'^ + 8M^e“^(l — T)aE2,a+i[du] 

Ao E £o -\- <A')E2Ad'^\ E E2^a+i[du\ + 

+ (^0 E £q + a^£o) (F2,a[c^w] + ||c?0|!^2,a+a) 


(5.49) 


1 


drY^WdO, 

c 

\\dO\ 


2 

C,fl 

< 11.7^112 


d ||2 
\ ^2 


+ 4M^e-\l - T)a'^\\dO: 


4\\2 

1^2 


'H" 


■ E 2 ,a+i[du] + £o{E 2 Adu] + \\dO\\ 




for all T G (0, T). 


Assuming Proposition |5.5] we prove the above claim (5.47). After summing (5.48),(5.49), 
we absorb into the LHS the critical terms 


E2,a+i[du], \\dO\\ 




which appear in the RHS of the above inequalities. This is done by picking the 


parameter a sufficiently large (but finite). The contraction estimate (5.47) then 
follows from Gronwall’s inequality for T > 0 suitably small. 


Proof of Proposition \5.5[ The proof follows exactly the lines of the proof of Propo- 
siti on|5.2[ The only notable difference lies in the estimation of the analogous term 
to (5.28), derived in (5.29)-(5.34). We sketch the argument in the present situation: 


Let J denote at most one spatial index, | J| < 1, either 1,2 or 3. The main term to 
be estimated is 


(5.50) 


(recall H < ^) 

r2 


g(J) 


RHS of ( |5H^ 


^do{dui,)ijj 


k.2q;—3 


< 


I do{du,,)ipj 


lL2 


I [RHS of 


^Q!—2 


|2 

Il2’ 
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where {du^) i^j := d^'^\du^)ij. Plugging in (5.43) and using the basic estimates in 
Proposition along with the assumption (5.46) we obtain 


(5.51) 


[RHS of ([^431 


pa —2 


|2 

Il2 


< 


(OPa a(du) + OTsdu + OPgdO) 


lL2 


{dO [Paati + Pau + Pa (O - /) + Pgao]) 


pa —2 


lL2 


I (or 3 a(dO)) ,,2 , „ [{o + o)dOd^{u,),,] ,,2 

I ^q-2 IIl2 

a(“')G(dM,dO) 


L2 


\\h 

I - 

(recall the asymptotics (|4.5|)) 


^ iia(‘^)a(du)||2 , iia(du)||2 , „ du ,,2 , ^ tii^(^^)ii2 

- II ^a-i + ll^^lli^ + II;:;i^IIl^ +^o(||^;^iIl2 


+ £0 


I 7* 2 

,a(‘^)(dO),,, ..dO. 


pa+ „ 


2 


-lli^ + ll—lli~(ll«ll?^=.= + llo-/|| 


^ 2 , 0 +f, 


.d^-^^didO) 




15(dO), 


lli^ + ll^lli= + ^oii^lli2 


I y-a+l 11-^ II ll-G II 

[(O + G)dOa 2 (u,)y] 2 „ d^J^G{du, dO ), 


^ -^2,a+l [du] + £oE 2 ,a[du] + ||dO|| 


L2 ^ 

2 


2 

L2 


+ £o\\dO\\ 




du 

rCk + 1 


\\h) 


+ £olMO|! 


I a^-^) [(o + o)doa2(u,),,] 2 „ a(^)G(du, do) ,,2 




lL2 


11^ 


We proceed to the problematic term (O + 0)dOd'^{u^)ij which can be controlled 
only in H^: 


(5.52) 


aG)[(o + 5)doa2(u,),,] 2 


<lia(")(o + o)||i.||do|L 


r.a-2 IIl2 

2 ||2 , iiTi , 7%||2 II I 


S^lli. + l|o + o||i. 


j.a-2 iir.2 ' II-.- I —iin^ii ^f-i Hi'* 


7* 2 


^11^4 + 110 + 0||i=.||- 




II2 


a —2 


lli= 


(employing the L'^ estimate (5.4)) 


< 


i£i + a^£^)\\dO\\ 




Finally, plugging in the nonlinearity (5.44), we have the bound 

ia(‘^)G(du,dO)||2 


(5.53) 


lL2 
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< 


(r 3du{u + u + + uu)) 


2 9(‘^)(Ou9(dM))||2 

o 


IL2 


L2 


I {Odudu + dOudu + Od{0)d{du)) ,,2 


pa— 2 


lL2 


5(-^) (^Od{dO)du + dOd{0)du) ,,2 


2 


lL2 


< (fo +£2)i?2.a[rfR] + (fo + i)IMu||i=ll^^lli. + 


{So + £S)\\dO\\ 


d{dO)d'^u 2 




ra — 2 


L2 + ||dO||_f^oo£o 


<{Eo+£o){E2,am + \\dO\\ 

, ii'9(c^O)„2 ,1 d^Tl „2 

+ II °_i IIl'iII “_i lli^ 

7* 2 7* 2 






, 5^0 


(by the estimate (|5.4| 


< 


(£"0 +£-0 0!^^o){E2^a[du] + ||dO|| 




□ 


6. The constraint equations in a singular background of unbounded 

MEAN CURVATURE 


In this section we prove Theorem |1.2[ our main stability result for the constraint 
equations The proof is based on mapping properties of the constraint map 

HO] 0 H and the implicit function theorem. Although similar results have been 
achieved in the smooth case and some rough backgrounds (see [7] for a general 
exposition), to our knowledge, the singular Schwarzschild background ((j^ eludes 
the standard references in the literature. 


It would be interesting to obtain perturbations (g, K) oi the Schwarzschild initial 
data set {^g,^K) on the singular hypersurface Eq, which satisfy the assumptions 
of Theorem o in a more constructive way. A general such construction would 
make use of the conformal method [3] that we do not employ here. However, one 
of the main obstructions to overcome in this approach is that the mean curvature 
tr-gK of the perturbation is unbounded. In fact, one can check ((|^ that 

tr^AT ^ LP(So), 

The results in the literature of the constraints using the conformal method were 
restricted in the past to the constant mean curvature (CMC) or ‘near CMC’ regime 
[3]- Recently, there has been a number of advances to the case of large mean 
curvature, ‘far from CMC ’, [nmani]. Yet these results contain certain regularity 
assumptions which in particular imply that the mean curvature of the initial data 
set is in L°° and therefore do not apply to our case. 


We consider below various tensors living in the weighted spaces iJ®’“(Eo) we have 
introduced (4.9). It will always be implied that this is understood by assuming 
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their components, evaluated w.r.t. the Schwarzschild frame (3.10), lie in the same 
spaces. 


6.1. The constraint map; Linearization and stability. Let 

(6.1) X ^ x (Sg), 

■^idg,j,dK,j) := (R - - VRr^X) 

denote the constraint map of the perturbed initial data 

( 6 . 2 ) dg,j := 5 ^ - ^g^j dK,j := K,j - ^Kij, i, j = 1 , 2 ,3, 

on Eo := (-e,e)^ x r(T = 0 )S^ 

Lemma 6.1. 'I' is well-defined, bounded and (Frechet). 


Proof. First, we express in terms of differences of perturbed and Schwarzs child 
components. Since the Schwarzschild pair {^g,^K) satisfies the constraints (1.5), 
we write schematically: 

(6.3) ^{dg, dK) = {d{dA) + {dd){^ A) + ^ AdA + {dAf + ^KdK + {dKf, 

d{dK) + {dd){^K) + ^KdA + ^ AdK + dAdK), 

where dA := A — ^A. The boundedness of 'I' now follows in part from Lemma 


using the asymptotics (3.15),(3.16), and by applying the L°° bound (5.3) on the 


4.5 


quadratic terms (dK)"^, (dAy, dAdK. Realizing that dA is in the variable dg (by 
expressing dA in terms of dg), we conclude that 4' is also in the above weighted 
spaces (taking into account the previous asymptotics and properties used). □ 


Remark 6.2. In view of the asymptotics (3.6) of the radius function r, at r = 0, 
we may replace the norm (4.9), at our convenience, with the equivalent norm 


(6.4) 


equiv 




|J| = fc<s • 


rp2a — 3{k — l) 


g-i 


where '®V is the induced Schwarzschild connection on (Sg,'^ g); we replace with 




in (4.9) as well, since the weights in the norm can exactly tolerate the most 


singular coefficients (3.15) of any additional generated terms. 


Our strategy is to fix the boundary components (close to Schwarzschild) and then 
solve for the variables in the interior via a perturbation. For this purpose we 
introduce the functions 

( 6 . 5 ) dg^°'> € r 4>“+5 n ny dK^°'> S H^’°‘ n dg^^'> G G 

dg = dg‘'°^ + dg^'^'f , dK = dK^°'> + , 

where by Hq we denote the set of functions in id^(Eo) having zero trace on the 
boundary of Eoj^ We fix dg^^\ dK^^^ and study under variations of the variables 


^^The boundary OSq has two components diffeomorphic to the sphere. 
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dg^^\dK^^\ Whence, we can view as a map 

^ : ((i74.a+i nH^) X X (ij4.“+i x ^ £r2.a-f ^ 


( 6 . 6 ) 


'i>{dg^°\dK''°\dg^^\dK^'^'>) := «'(d5,dii:). 


Let 

(6.7) := 7?(dg(o).dif(o))«'(o,0.0.0) : x ^ x 

be the Frechet derivative of w.r.t. the variables {dg^^\ about zero, that is, 

the Schwarzschild initial data; dg^^^ = dK^^^ = dg^^^ = dK^^^ = 0. The following 
theorem is our stability result for the constraint equations. 


Theorem 6.3. Let a > 0 &e sufficiently large, eonsistent with Theorems \4. ?! \4^ 
Then Z?d> is surjective, and the level set 

-1 


( 6 . 8 ) 




({ 0 }) 


(•.•.dg(i).dif(i))y 

is non-empty, for every pair {dg^^\dK^^^) in a sufficiently small ball of H 4.“+2 X 

H3,a 

Remark 6.4. The preceding theorem implies the existence of non-spherically sym¬ 
metric perturbations (g, iLy )[^of the Schwarzschild initial data set {^g,^Kij) on 
Eq, compatible with the constraint equations (1.51, which lie inside the weighted 
spaces i74’“+2(Eg) x satisfying the assumptions of the local-in time- 

existence Theorems 14.7114.91 


Remark 6.5. We can construct perturbed initial data sets {g, K) verifying the 
constraints (1.5) and the assumptions of Theorems 4.7 4.9 which are identical to 
the Schwarzschild initial data set {^g, ^K) on Eg H {a: < 0} and not spherically 
symmetric on the other ‘half’ Eg n {a: > 0} (and vice versa) by applying the 
arguments in this section only to the {x > 0} ‘piece’ of the initial data. 

Unfortunately, the formula of Hdt is quite complex [20] |8] and it is not easy to prove 
directly the surjectivity of H'l'. We do not present the formula of H'L here, since 
we are not going to use it. Instead, we follow an indirect argument of |5], which we 
adapt in our context. We consider the restriction of (6.1) in the subdomain 

K,,=ip^{^K + Cx^g),,, 


(6.9) 


9 = 


- 


where {Cx°9)ij '■= -\- ^VjXi — ^gij^V X^. Then the constraint map (|6.l|) 

becomes (0) 


^ : (ij4.«-Hi n Rgi) X (ij4.a-Hi n Rgi) ^ ij2.«-i x 


(6.10) 4'((/5-1,X) = (v? ^[-8Asg(p-\-(p^R-yy\^K-\-Cx^g\l- 


+ Cx^g)^j + 4ip-\^K + Cx^gffi^Xjip 
-2y,-^^VMt^s/K + Cx^g)]) 


^^This is easily seen by considering non-spherically symmetric 
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Further, the linearization of ^ at (0,0) reads 

Y) = {- SAs-gr, + 2{tis/K)^v - 2K^^iCY^g^ + trs/ 

( 6 . 11 ) ^V\CY^9h + - 2tTs/K^W,r] - 2^V^^K,,ri) 

Recall (3.151, (3.6) to find the leading asymptotics, as cc —)■ 0+, r = 0, of the 


singular Schwarzschild components of D'^{ri,Y): 


( 6 . 12 ) 


’if,,I, |trs/if| 


< 


X2 




1 


X2 


(|J| + 1) 


A key ingredient in the proof of Theorem |6.3| is the following proposition. 

Proposition 6.6. The operator 

(6.13) oil : iih«+i n 

is Fredholm. In particular, it has finite dim kernel and cokernel. 


Assuming Proposition |6.6[ we proceed to the proof of Theorem |6.3| 


Proof of Theorem \ 6. 5| Evidently, rangeH'I' D rangeUT and hence by Proposition 
6.6 HT has finite dimensional cokernel. It follows that rangeU^' is closed. 


Since rangeil^' is closed, rangeil^' = (kernelil^'*)-*-; D'^* being the adjoint of . 
Therefore, in order to prove surjectivity, it suffices to show that kernelH^'* = {0}. 
We now recall a lemma from [20]: 

Lemma 6.7. The kernel of D’^* is in one-to-one correspondence with the set 
of Killing vector fields of the Einsteinian vacuum development of{g,K), i.e., the 


Schwarzschild region (3.1) foliated by {S.r}TG[o,T]• 


Picking a large enough, we can guarantee that none of the Schwarzschild Killing vec¬ 
tor fields {dt, d^, sin (fdeP cot 9 cos (j) d^, cos (j)dg — cot 9 sin (f d^) lie in 2 (Eq), 

the domain of H'k*. Therefore, for this a the condition kernelH'I'* = 0 is verified 


and hence is surjective. The second part of Theorem 6.3 follows by the implicit 
function theorem. □ 


6.2. Proof of Proposition 6.6, We realize the following plan (which is an adap¬ 
tation of a usual argument from the non-singular case): Given cr > 0 we define the 


bounded operator (see Lemma 6.1) 
(6.14) 


5 := I?# + — (/„, -Iy) : n H^ 


H 


2,a-% 


where Iri,lY are the identity maps for the variables r],Y respectively. Then we 
show: 

Proposition 6.8. S is an isomorphism for a > 0 sufficiently large. 


Using the preceding proposition, we obtain the Fredholm property claimed in 
Proposition |6.6| as follows: 

Since S is an isomorphism (for large a > 0), it has a bounded inverse. In fact. 


( 6 . 15 ) 


5-1 : i/2,a-i ^ i/4.a+i ^ 
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is a compact operator by a weighted version of Rellich’s theorem. By definition 


(6.14) the operators 


(6.16) nil oS-'^ - I = -^{In,-lY)oS-^ ^ 

S-^oD^-I = - S'"^ o -Iy) : n i/ha+f p 


are also compact, since the RHSs of (6.16) consist of compositions of bounded with 


compact operators. Thus, HT is Fredholm. 

A note is in order here. The key ingredient in the previous proof is the invertibility 


of the operator S (6.14). As we shall see below we are able to prove Proposi¬ 


tion 6.8 thanks to the precise analogy in the leading blow up orders (6.12) of the 


Schwarzschild components appearing in the expression of the linearized map 


( 6 . 11 ). 


We divide the proof of Proposition 6.8 in two steps. Recall the formula (6.11). 


Proposition 6.9. For every h, Zi € H^’°‘ ^, i = 1,2,3, the system 


(6.17) -SAsgT] + 2{trsg^K^ri - 2^K g)i] + trsg^K{^\I Y)i H—= h 


SvyJ I 


V{CY^g),g + 4*)a/^V,?7 - 2trs/K^V,g - 2^V^K,gg - = 


■)S\-r3 S 7 


has a unique weak solution g,Y G Hq 


l,a+f 


, provided a > 0 is sufficiently large. 


Proof. Define the bilinear form 


(6.18) B : X H, 


l,a+% 

0 X ^0 


S[(r7,P),(?,X)] = (LHS of ([^,( 


A 


' ^2q:- 1-3 ^ ^Sck-j-S ' 




We want to apply Lax-Milgram. For that we must establish the two estimates: 


(6.19) \B[ig,Y),i^,X)]\< 
and 


p,„.3+||y||2 )(M_3+||Af 


( 6 . 20 ) 


B[{v,Y),ffi,Y)]>Ml,^^^, 




The bo und (|6.19 ) follows easily by Cauchy-Schwartz, taking into account the asymp¬ 
totics (16.12^ of the coefficients of (16.171). Hence, it suffices to show (16.201). We 
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proceed by integrating by parts, employing the asymptotics ( 6.12| ): 

rj Y 


S[(77,y),(r?,r)] = (LHS of ([^,( 


/Eo ^ 


2q;-1-3 


— Ca 




- c 




/So X 


201+3+ o 


d/is 


/So X 
|Syy|2 

,y»2o+3 


2o+3+ 7 


dfiSq + a 


77" 


/so a/" 


))l^> 


rd/iS 


d^iSg — 






- Ca 


- C 


++Vani. 

+>+,771 


2o+3 

y2 


d/is 


/So a; 
'S 


/So a/ 


2o+3+ 7 


d^Sg — C 


So ^ 


2o+6 


d/i6 


/So ^ 


2o+5 


(utilizing the estimate |'®Rfci| < x ^ and applying C-S) 


’V77P 1 

-d/rsg + - 


’vrp 

1^20+3 


d/isg + (cr — Ca^) 


+ 


p2o+6 


d/is^ > 


d/is 


y2 


/So 2; 


2a+6 


d/is 


Taking cr large enough, a > a^, we arrive at (6.20). Thus, the Lax-Milgram theorem 
can be applied for the system (6.17), from which the conclusion of the present 
proposition follows. □ 


Remark 6.10. Away from x = 0 the coefficients of (6.17) are smooth and bounded 
(locally). By standard interior elliptic regularity theory we derive that the weak 
solution in Proposition 6.9 r],Y lies in Hf^g{Yo \ {a; = 0}). 


Next, we show that the solution of (6.17) is in fact ri,Y G This implies 


that S (6.14) is one to one and onto, proving Proposition 6.8 

usiLiuii ± iit^ wcaK suiUL'iuii ij, 1 t rl ’ 

I ., q 

sition 


Prop osit ion 6.11. The weak solution r],Y G Hq 
satisfies r],Y £ 


6.9 


of [6.11) furnished by Propo- 


Proof. rj^Y G iJ^’“+2; Multiplying (6.17) with (-^+, —^5+), commuting covari¬ 
ant derivatives and integrating by parts twice we arrive at 

f |^V^V277|+ ^ f |Va2 + V277l , f 


/So 

s 


2a -dTSg - Ca 


S\ 


y,2Q;+ ^ 


-dfis 


X 


2cx. 


+ 22 / 7 1 

t^2q;+^ 


< 


g.2a+2, 

So 2/ ^ 

I+V2277I 


dfiS g , 


I^v/V2y+ 

rri2ClC 


djasg — Ca 


+a2y+V2yz| 

p2q;+ i 


I^V 2 y*P 


JHq X 

Riem- — - d^is-r 


X 


2ol 


+ ^ '—fi£^dp.Sg 

•/So 2/ 

f / '=’Riem- ^ —d/rsg 


- C 


’Va + V22T.| 


X 


2a-|-| 


dpLSg — C 


| + V22T^| 


^2q:+2 


dflSg < 




dfJLS 
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Using the bound '^Riern < ^ and applying Cauchy’s inequality we deduce for a 
large (o' > a^): 


|®V^V2ry|2 


d/is 


< 


'H 




riiL 


/i2 






d/is 


I 

-E 




d^s 


X 


< + 00 , 


since h^Zi € 2. Likewise for the d'i (i.e., the other rotational d irecti on). The 

■^Vii derivative of the functions 77,!) is estimated separately: From (6.17) we have 

""ViiDp , , f I^VnryP 


E 


X 


2ol 


-dfxs 


x^ 


^E 

a— 2,3 


iSyjS 


r.2a 


■dfiSg + 


’VaaT?P 

/-p2oi 


dfis 






+ cr\\v\\]ji ,^+3 + W^Wh +'^\\^\\h < + 00 . 

i 

Thus, 77,1) S . The full estimate is obtained by differentiating the 


system (6.17) in the rotational directions 92,i93, applying the same argument with 
the analogous weights and then/ differentiating (6.17) in the di direction. □ 


Appendix A. Changing frames freedom; Propagating identities; 

Retrieving the EVE from the reduced equations 

Given a spacetime g) and an orthonormal frame {ei/p, one may change to a 

Lorentz gauge frame {ei}o by solving the following semilinear system of equations, 
which is derived by taking the divergence of (2.9): 

□g(Oi) = {diYAx)a‘'o\ + AdO + AdO + 0^,diy{Ax)k 


(by (2.15) for A) 


(from (2.8)) 


= A'^0 + AdO + AdO + OdivA 
= O^A^ + 0^{dOf + AO^dO + AO'^dO 
+ 0{dOf + AdO + OdivAjv, 

where the terms without indices in the RHS stand for an algebraic expression of a 
finite number terms of the depicted type. 

Lemma A.l. If the above system (which we write schematically as) 

(A.l) ^g{0\) = O^A^ + 0^{dOf + AO^dO + AO'^dO 

+ 0{dOf + AdO + OdivA. 

is well-posed in a certain solution space, then there exists a unique orthonormal 
frame 

(A.2) g, = 0/efa 

with 0\ lying in that particular space, which is identical to {ei/p on the initial 
hypersurface Ep, verifies the Lorentz gauge condition (2.15) and such that the con¬ 


nection coefficients (Ap)ij := g(Vggei,ej), i < j, are equal to a priori assigned 
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functions on Sq; within the corresponding space of one order of regularity less than 

Ol 


Proof. It suffices to show that the initial data for (A.l) is uniquely determined by 
the assertions. We set 


(A.3) 0\{t = 0) := 

Let 

(A.4) eo(0'’)(T = 0) = eo(0'‘)(r = 0) =: 

Then the transition formula ( |2.8[ ) for X = oq reads 
(A.5) (Ao)jj(r = 0) = (Ao)jj(r = 0) + 


(i.e., e* = Ci on Sq). 


Thus, the components (Ao)i_; can be freely prescribed initially by choosing h\ in 
(A.41 accordingly. □ 


A.l. Proof of proposition |2.ll We will leave the reader to fill in the details for 
the fact that the solution {A^)ij,Of of (2.161,(2.19) corresponds to a spa cetime 
(/). This is a consequence of the necessary initial assumption (2.20). One 


such immediate consequence follows from (2.19) for z = 0: 


(A.6) 


5o(Oo) = -oln 


0-^ [Ob] > 


Og(r = 0) - /o“ = 0, 


which implies Oq = /o° and hence eo = do everywhere, since = 0. The set 


of functions 0“ defines the orthonormal frame {b^Iq in through (2.18) and 


hence completely determines the metric g. What remains to be verified is that the 
connection coefficients of {cijo are indeed the {A,y)ij's of the given solution. In 
other words, we have to show that the connection D induced by the solution set 


(A.7) 


• - (.Ai/')i C/j;, 


is the Levi-Civita connection V of the metric g. Formally, one cannot take this 


for granted. It has to be retrieved from the equations (2.16),(2.19) and the initial 


assumption (2.20). For example, the compatibility of D with respect to g is encoded 


in the skew-symmetry of the {AX}ij's 


(A.8) 


D{g) = 0, 


iff 


-j- — 0, 


which also has to be verified, since it is a priori valid only initially (2.20). The way 


to do this is by deriving the following new system of equations from (2.16) for the 
symmetric sums: 

ek{{Aff)ij -\- -j- \{Ak)ij -\- {Ak)ji\) 

(A.9) -I- By [A[{A)ij -\- {A)ji ]) -I- + {Ak)ij \), 

where we have assumed that the sum {A‘^)ij+{A'^)ji corresponding to the ter m JC 
in the gauge condition (2.15) can be expressed as A[(A)y -j- (A)ji]. Since (A.9) 
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has zero initial d ata (|2^2^ , the symmetric sums are zero everywhere and hence the 
skew-symmetry ( A.8[) propagates 



vanishes and we are left with 
(A. 10) □g(divAx-A2)^^. 

= (A‘')j°e^(divAx - A'^) + {A'^)j’'e^{divAx - A'^). 


The Lorentz gauge condition is valid initially (2.21). If the cq derivative of (divA^ — 
is zero as well on Eq, then the Lorentz gauge is valid in all of = 

E X [0,T]p^ This is in fact implied by (2.17), putting i/ = 0 we have 


(ATI) 


eo(divAx - A"^) = V^Rq* - ViRoj = 0 


on En 


by virtue of the vanishing of Ra6(T = 0) (2.21) and the (twice contracted) second 
Bianchi identity, V“Rab = ^R, to replace if necessary a transversal derivative with 
tangential ones to Eq. 


On the other hand, taking the V® divergence of (2.17) and commuting derivatives 
we obtain 

1 . 


lEgRjyj — V®VjR^i — — VjVj/R-l-R’ 


J V 


R^ 


R® R 


(A.12) 


= R%- .,Rc 


R/Ri. 


where we employed again the twice contracted second Bianchi identity and the fact 
that the scalar curvature R vanishes everywhere: [contracting {vj} in (2.17)] 

(A.13) 


0 = ViR- ^V,R= ^V,R, 


RL = 0. 

I 2^0 


Now that we know the Lorentz gauge is valid, the identities (2.21) and (2.17) i = 0 
imply 

(A. 14) R^^- = 0, VoRi/j = VjR^O) on Eg. 

Utilizing the second Binachi identity V“Ro& = ^R = 0 once more, we conclude that 
VqRi/j vanishes and hence the initi al dat a set of (A.12) is the trivial one. Thus, the 
initial condition = 0) = 0 (2.21) propagates and the spacetime 

obtained from the solution of ( 2.16| verifies the EVE (1.1). □ 

Remark A. 2. Given the frame |ei}o initially on Eg, and once the components 
(Ag)ij(T = 0) have been chosen rd then the initial data set of (2.13) is fixed by 


20r 


This follows by a basic a priori energy estimate for linear systems like (A.9i, which in the 


singu lar Schwarzschild background is derived in §5.3| for the more involved quasilinear system 

(O. 

^^Note however that the term eo(divAx — of second order in A and hence not at the 

level of initial data for ( |2.16| l, which we are allowed to prescribe. If zero initially, this should be 
a consequence of the geometric nature of the equations. 

^^The (Ao)ij’s are not fixed by the Lorentz gauge condition; cf. Lemma . 


A.l 


They correspond 


to the ^0 derivative of the frame components e^, which we can freely assign initially. 
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condition (2.21), i.e., the EVE and Lorentz gauge on Eq- Indeed, the components 
= 0)j = 1)2,3, are determined uniquely by the orthonormal frame 

{cilf on (So,g). The (A^)oj(t = 0)’s correspond to the components of second 
fundamental form Kij of Eq, which is given by the solution to the constraints (1.51, 
included in (2.21 ). Moreover, the expression of (2.21 ) in terms of for z = 1, 2,3, 
reads (schematically) 

(A.15) eo(Aiy)oi = ei/(A) H- 


^o{Ao)ij — ^ ^ e^(A^) 






-A" 


on Sn 


Hence, the LHS functions are expressed in terms of already determined components. 
Finally, the rest components eo(Aj^)ij(T = 0), z,j = 1,2,3, are fixed by the 
algebraic property of the Riemann tensor 


Hoi/ij — ^ijOiy 

^o(Ai/)ij Ci/i^Ao^fj ([Ayj,, A^])2j (A|'^)y| (A/j)^^ 

^i(Aj)oi>' ^ji^Ai^Qjj ([Aj,Aj])o^' (Aj'2)jj (A;^)q^. 


or 

(A.16) eo(A^.)jj — C(y(Ao)ij Ci^Aj^oj ej(A^)oi>' A , on Sq, 


since all the terms in the RHS have been accounted for. Notice that the dehni- 
tion of Riemann curvature was implicitly used in deriving (A.121 upon commuting 
covariant derivatives. 
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